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1. INTRODUCTION 

In a wireless network, relay station advantages can be evaluated through two parameters: performance and 
cost. From the performance side, relay stations can be utilized to extend the achivable data rate within the same 
cell or alternatively, they can be used to extend the coverage area [1]. When a relay station is installed to extend 
the coverage area, both, relay station and base station, use the same frequency at the same time which increases 
the spectrum reuse. Thus, the use of relay station improves the overall system throughput. By installing more base 
stations instead of relay stations, the same or in fact better performance can be achieved. However, installing base 
stations is much more expensive than installing relay stations. 

Relay technology can be used in rural scenarios to extend the coverage [2]. It can be used in the case of 
earthquake or disasters where deploying a fixed line backhaul link for a base station is difficult. In the last decade, 
cooperative diversity strategies using randomly distributed relay nodes between the communicating terminals have 
been extensively studied as their improvements in performance do not require additional power or frequency spec- 
trum [3-14]. The main objective of this work is to propose efficient relaying techniques to increase the sum rate 
and reduce the symbol error rate (SER) with low computational complexity. 


2. SYSTEM MODEL 

Let us consider a half duplex system which consists of M single-antenna mobile stations (MSs) commu- 
nicating with another M single-antenna mobile stations via a relay station (RS) having N ( N > 2M) antennas 
as shown in Fig. 1. There is no direct link between mobile stations and their communication partners. Relay 
station uses either the decode-and-forward (DF) or the amplify-and-forward (AF) protocol depending on the used 
technique. The noise at the relay station and at MS nodes is assumed to be modeled as independent, zero-mean, 
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complex Gaussian random variable with variance Tae and o”, respectively. Let us assume that the uplink and 
downlink channels are reciprocal and frequency flat fading. Further, the channels are assumed to remain constant 
during each transmission cycle. The maximum transmission power at the relay station and at the it” MS, i.e., 
MS;, is given by Prs and P;, respectively. Furthermore, it is assumed that the channel state information (CSI) is 


available at the relay and mobile stations. 
© 


RS 


N antennas 


, 
hu 


Figure 1. System model. 


The combined multiple access channel H € C’ <?™ from all MSs to RS is given by 
H=(h, by .. hy bY bh... hy] (1) 


where h; and h/, i = 1,..., M are column vectors representing channel from MS; to RS and from its corresponding 
partner, i.e., MS‘, to RS, respectively. Similarly, broadcast channel from relay station to all MS nodes is given by 
HË c C?MXN, For the given system model, in the first time slot, all users transmit their data to the relay station. 


The signal received at relay station yp € CN *! is given by 
M M 
Yr = X his; +X his, +npg (2) 
i=1 i=1 


where s; and s’, are the signals transmitted from the MS; and MS’, nodes, respectively, to the relay station and ng 
is the noise vector at the relay station in the first time slot. This signal needs to be processed in order to mitigate 
interference and noise. Let G € CYX be the processing matrix at the relay station. In case of using the AF 
protocol, relay processing matrix is represented by a single matrix, i.e., G [15]. Whereas, in case of using the DF 
protocol, relay processing matrix is represented by a multiplication of three matrices Gm, Gp, and W, such that 


G = GmWG, (3) 


where Gm is used to remove the effect of the interference occurring in the multiple access phase, i.e. during the 
first time slot, permutation matrix W is then used to rearrange the resulting vector in a proper order before sending 
it, and G, is used to remove the effect of the interference occurring in the broadcast phase, i.e., during the second 
time slot. After processing the received signal y p defined in (2) at the relay by using relay processing matrix G, the 
signal vector Xp is obtained which is then transmitted to all mobile stations in the second time slot. 


3. TECHNIQUES TO MITIGATE INTERFERENCE 

This section presents proposed relaying techniques to mitigate interference. The orthogonal projection 
technique enjoys low computational complexity, however it suffers from low performance in terms of SER [15]. 
The other technique is based on maximum likelihood (ML) strategy to detect the symbol vector and then uses 
minimum mean square error (MMSE) strategy to broadcast the resulting vector. This technique enjoys optimal 
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performance in terms of SER, however it suffers from high decoding complexity due to the use of ML detector at 
the relay. More techniques are also proposed in this section in order to reduce the overall computational complexity 
and improve the overall system performance using network coding and randomization strategy as explained in the 
next subsections. 


3.1. Multi-antenna technique based on ML and MMSE strategy 


In the technique proposed in [15], a zero-forcing strategy is used to reduce the effect of the interference 
at the cost of noise enhancement [16]. In order to improve the previous technique, other strategies can be used to 
reduce the effects of interference. Note that both, the interference in the multiple access phase and in the broadcast 
phase, need to be mitigated. In this technique, ML detector explained in [17] is used to detect the received signals 
at the relay during the first time slot. To mitigate the interference occurring in the broadcast phase, MMSE strategy 
is applied. 


3.1.1. ML detector 


As explained in Sec. 2., we are considering a multi-user system where all users are transmitting their 
signals at the same time to the relay station and the received signal vector at the relay station yp is given by (2). 
The key idea of ML detector is to find the joint error u for each possible combination of the transmit symbols, such 
as 


N M 
u=} lya(n) - do lhiln)si +h;(n)s;)|? (4) 


where yr(n) is the nth element of the vector yr. After calculating the value of u for each possible combination 
of the transmit symbols from all mobile stations, the detected transmit symbol vector § € C?M X1 which gives 
minimum value of u is obtained. ML provides an optimal solution to the detection problem, however it suffers 
from extremely high decoding complexity due to the exhaustive search over all possible combinations of symbols 
where its decoding complexity increases exponentially with the increase of the constellation size and the number 
of transmitted symbols. 


3.1.2. Permutation matrix 


After detecting the symbols optimally using ML detector, the symbols need to be arranged in a vector in a 
proper sequence in order to receive them correctly at the destination node. The permutation matrix W € C2” *?™ 


used in (3) is given by 
w= Ea Im ) (5) 
Iv Ouxm 


where Oy7x x4 and Imxm denote an M x M matrix which contains zeros in all its entries and an M x M identity 
matrix, respectively. The detected symbol vector § is multiplied by the permutation matrix W, such that 


t= Ws. (6) 


3.1.3. MMSE strategy 


To mitigate the interference in the broadcast phase, MMSE filter is used. This filter is represented by the 
matrix G, used in (3). After normalization, it is given by 


Gs = aG, (7) 
where Gp is the MMSE filter, given by 
Gs = H(H”H + 1/7In)', (8) 
and a is the factor to fulfill power constraint at the relay station, given by 


Prs 


a= > (9) 


Vič? 
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The transmitted signal vector from relay station can then be obtained as 
Xp = Got. (10) 


The signal given by (10) is then received by all mobile stations where the signal received at the it” mobile station 
MS; is given by 

yi = h” xg + ni. (11) 
In this technique, MS; needs only its CSI, i.e., h;. The achievable sum rate of this system can be calculated by 


using 
2x M 


i=1 
and the receive signal to interference and noise ratio (SINR) at the MS; in this case is given by 
[h Goi? 
H 
o? + paren |h; Go; |? 


y= (13) 


where Gp; represents the it? column of the relay transmit filter G+. 


3.2. Multi-antenna technique based on MMSE strategy 


In this technique, MMSE filter is applied at the relay twice. First, as a receive filter at the relay during the 
first time slot. Second, as a transmit filter before sending the decoded data from the relay station in the second time 
slot. During the first time slot, the MMSE detector at the relay is given by 


Gm = (H”H 41/71) H”. (14) 
Using this detector, the estimated signal vector § is obtained as 
S = Gp. (15) 


This estimated vector is then hard decoded and used for further processing. Note that the estimated vector without 
hard decoding could also be used. The decoded symbols need to be rearranged to make sure that every node 
receives symbol from its corresponding partner. This can be achieved by using permutation matrix given by (5). 
The rearranged vector t is given by (6). Before transmitting this signal vector, transmit MMSE filter given by (7) 
is performed. Sec. 3.1.3. also explains the calculation of the power normalization factor a expressed in (9) for this 
filter. The final transmit vector xp is given by 

Xr =Gyr (16) 


where G is defined in (3). The achievable sum rate using this technique can be calculated using (12). 


3.3. Multi-antenna technique based on network coding 
3.3.1. Concept 


Network coding can be useful in combining signals to transmit them using less number of time slots [3, 
7-12]. In multi-antenna scenarios, network coding combines signals to reduce the number of transmitted symbols 
which can then be transmitted over less number of antennas. Otherwise, if the same number of antennas is used 
to transmit these combined symbols, a better performance can be achieved. In this technique, M-PSK modulation 
scheme is used. Making use of the fact that when two M-PSK symbols lying on the unit circle are multiplied, the 
resultant symbol lies on the same circle. Because of this property, E {ss } is preserved even after multiplication is 
performed. 


3.3.2. Implementation 


As explained earlier, the basic idea behind this technique is to use network coding at the relay to improve 
the performance in terms of sum rate as well as SER. The received signal vector at the relay station is given by (2). 
This received signal is the sum of all the signals coming from all the mobile stations. A receive filter is needed 
at the relay station to separate the signals and thus detect the correct signals from all the nodes. In this technique, 
MMSE filter given by (14) is applied. The estimated signal vector § is obtained using (15). The symbol vector $ 
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is then hard decoded based on the decision boundary of the used modulation scheme. Note that the signal vector $ 
can also be sent without hard decoding. In the next step, symbols belonging to the same communication pair are 
multiplied to obtain a new combined symbol. The new symbol tų, generated from the symbols of the kt” pair, is 
given by 

th = k Fh k=1,....,M. (17) 


To keep the symbol t generated for the kt” pair separate from the symbols of the other pairs, the orthogonal 
projection strategy is applied. The key idea of the orthogonal projection strategy is to find a precoding matrix at 
the relay that groups the signals from the same pair together and eliminates the inter-pair interference [15]. The 
constraint for finding such a precoding matrix is N > 2M — 1 which is fulfilled by our system model. In this 
strategy, MS; needs to know the relay precoding matrix in order to remove the effect of self-interference which 
means that the CSI is required at the mobile stations. The precoding matrix for the k*” communicating pair is given 
by 

P; = [hy...bg—ihyy1..-byh)...by_ h-hh, | (18) 


where P, € CX *(2“—2) is a submatrix of channel matrix H, given by (1). P;, is obtained by removing the kt” and 
(k + M)* columns from H. The orthogonal projection matrix Q, € C’*% is obtained from interference channel 
P% and is given by 

Q; = (In — P (PEP) IPE) k=1,...., M. (19) 


Q; is then multiplied by the received signal y p given by (2) to recover the signals of the kt” pair, such as 


Note that Qh; = 0 and Q,h; = 0(k 4 i). Let us define a matrix 1 € C*1 which contains ones in all its entries. 
The transmitted signal from the relay is given by 


M 


xr =a) Qa -1-t (21) 


k=1 


where a, used to normalize the transmit signal power at the relay station in order to fulfill the relay station power 
constraint, is given by 


(pea _ (22) 


VEz lOl? 
“th 


The received signal at the 7°” node MS; during the second time slot is given by 
Yi = h? XR +H Ni. 


In this technique, MS; needs only its CSI, i.e., h;. The achievable sum rate of the system can be calculated using 
(12) where SINR at the node MS;, is given by 


N H 
= [oper ah; Qi | 
a 2 


Vk (23) 


oO 


and Qx; is the j+” column of the relay precoding matrix Q, for the kt” pair nodes, i.e., MS; and MS’ x. 


3.4. Randomization techniques 
3.4.1. Concept 


In the techniques explained in Sec. 3.2. and Sec. 3.3., the MMSE strategy is performed at the relay to 
mitigate the interference occurring in the first time slot using (14). The techniques based on MMSE strategy 
with or without network coding are simple and enjoy a low decoding complexity, however, they suffer from low 
performance in terms of SER as compared to the techniques based on ML strategy. To improve the performance of 
the techniques based on MMSE strategy during the first time slot, let us search randomly for a symbol vector having 
same mean and variance as our estimated MMSE symbol vector defined in (15). In other words, the randomization 
strategy is used during the first time slot to find a better estimated symbol vector at the relay and thus to reduce the 
overall error rate. 
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The signal received at the relay y p is given by (2) and the MMSE filter G,, used for obtaining the estimated 
symbol vector § is given by (14). This estimated vector § is obtained by (15). After receiving yp and obtaining § 
using MMSE filter, a new random symbol vector §,,¢ is generated with mean and variance equal to the original 
estimated symbol vector $. Afterwards, two error values for the two estimated symbol vectors, § and new, are 
generated using 


e = |lyp —Hs\||’, (24) 
Yr — H8newll?- (25) 


If €new < e€ is true, this means that we have found a better estimated symbol vector. Therefore, the original 
estimated symbol vector § can then be discarded and the new one, i.e., ŝnew, can be stored instead. If the original 
estimated symbol vector S$ is better, then the newly obtained one, i.e., Snew is discarded. We again try to find an 
estimated symbol vector better than the obtained one in the last step. For that, we again follow the same procedure 
of finding a new random symbol vector, comparing it with the existing one, and then storing the best one. This 
process needs to be repeated for a pre-defined number of iterations to obtain a better estimated symbol vector. 
The randomization technique based on MMSE without using network coding is explained above and named as the 
randomization technique based on MMSE strategy. This strategy, i.e., randomization strategy, can be combined also 
with the techniques based on network coding explained in Sec. 3.3. in order to improve the estimated symbol vectors 
of each pair before combining them using (17). The latter technique is named as the randomization technique based 
on network coding. 


Enew 


4. RESULTS AND DISCUSSION 

This section presents the simulation settings and the results obtained for the techniques described in Sec. 3.. 
The performance of the proposed techniques and the one proposed in [15] are compared using two parameters: 
average achievable sum rate and symbol error rate. A Rayleigh flat fading reciprocal channel is assumed for uplink 
and downlink communication. Channel vectors are assumed to be independent and identical distributed (i.i.d) and 
remain constant during the whole transmission cycle. The whole CSI is assumed to be available at the relay station 
while mobile stations require their own CSI as explained in Sec. 3.. For all simulations, the number of mobile 
station pairs intending to communicate with each other is set to M = 2 and number of antennas at the relay station 
are set to N = 4. Transmit signal to noise ratio (SNR) is varied from 0 dB to 30 dB. The power at each mobile 
station is P; = 1, Vi. The power at the relay is assumed to be proportional to the number of antennas at relay, 
thus, Pr = N = 4. The noise power is assumed to be changing in accordance with the transmit SNR requirement. 


4.1. Symbol Error Rate 


To find the SER at each SNR value, 500000 symbols are transmitted from each MS. In the techniques 
based on randomization strategy, 30 iterations of randomization are used. Fig. 2 shows the performance of each 
technique discussed in Sec. 3. and the one proposed in [15] in terms of SER. In Fig. 2, the legend OP, MLtxMMSE, 
MMSEtxMMSE, MMSESigMulOP, MMSErandSigMulOP, and MMSErandtxMMSE denote the technique pro- 
posed in [15], the technique based on ML strategy explained in Sec. 3.1., the technique based on MMSE strategy 
explained in Sec. 3.2., the technique based on network coding explained in Sec. 3.3., the randomization technique 
based on network coding explained in Sec. 3.4., and the randomization technique based on MMSE strategy ex- 
plained in Sec. 3.4., respectively. It is clearly visible that the proposed technique based on the optimal ML detector, 
denoted by MLtxMMSE, enjoys the best performance and outperforms the other techniques, however, as explained 
in Sec. 3.1., it suffers from extremely high decoding complexity. On the other hand, the proposed randomization 
technique based on MMSE, denoted by MMSErandtxMMSE, enjoys low decoding complexity with almost the 
same performance as compared to the one based on the optimal ML detector. Moreover, the randomization strategy 
can surely improve the performance drastically in a comparatively less complex way. 


4.2. Average achievable sum rate 


Fig. 3 shows the performance of the proposed techniques and the technique proposed in [15] in terms 
of achievable sum rate where the legend OP, MMSE, and SigMulOP denote the technique proposed in [15], the 
technique based on MMSE strategy which is explained in Sec. 3.1. and Sec. 3.2., and the technique based on 
network coding which is explained in Sec. 3.3., respectively. As explained in Sec. 3.3., the technique based on 
network coding reduces the number of transmitted symbols by combining them, therefore, Fig. 3 shows a significant 
gain in the sum rate achieved in the case of the technique based on network coding as compared to the other 
techniques. 
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Figure 3. Average achievable sum rate vs SNR (dB). 
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5. CONCLUSION 

In this work, the performance in terms of SER and achievable data rate of a variety of newly multi-antenna 
techniques used for multi-user two-way wireless relay networks are proposed and studied. The multi-antenna tech- 
niques based on orthogonal projection enjoys low computational complexity. However, the performance penalty 
associated with them is high. Moreover, the techniques based on ML strategy enjoys high performance in terms 
of SER, however, they suffer from very high decoding complexity. On the other hand, the proposed randomiza- 
tion technique based on MMSE strategy enjoys low decoding complexity with almost the same performance as 
compared to the technique based on ML strategy. 
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